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Abstract: By the Pauli exclusion principle no quantum state can be occupied 
by more than one electron. One can put it as a constraint on the electron den- 
sity matrix that bounds its eigenvalues by 1. Shortly after its discovery the Pauli 
principle has been replaced by skew symmetry of a multi-electron wave function. 
In this paper we solve a longstanding problem about the impact of this replace- 
ment on the electron density matrix, that goes far beyond the original Pauli 
principle. 
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1. Introduction 

The Pauli exclusion principle, discovered in 1925, claims that no quantum state 
can be occupied by more than one electron. In terms of the electron density 
matrbflp this amounts to the inequality {tp\p\%p) < 1, that bounds its eigenvalues 
by one. The following year Heisenberg and Dirac replaced the Pauli principle by 
skew symmetry of a multi-electron wave function 11, Ch. 4]. 



1 There is no agreement on a proper normalization of the one-electron matrix. To avoid 
a confusion we call it electron density for Dirac's normalization to the number of particles 
Trp = N, and reserve the term reduced state for the probability normalization Trp = 1. 
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The subject of this study is the impact of this replacement on the electron 
density matrix. The latter determines the light scattering and therefore quite 
literally represents a visible state of the electron system. The impact goes far 
beyond the original Pauli principle. As an example, consider three electron sys- 
tem /V'TLq with one-electron space TLq of dimension 6. Then the spectrum A of 
the electron density matrix, arranged in non-increasing order, is bounded by the 
following (in)cqualities discovered by Borland and Dennis |3J 

Ai + A 6 = A 2 + A 5 = A 3 + A 4 = 1, A 4 < A 5 + A 6 . (1) 

The authors established the sufficiency of these constraints and referred for a 
complete proof to M.B. Ruskai and R.L. KingsleyO It worth reading their com- 
ment: 

We have no apology for consideration of such a special case. The general 
N -representability problem is so difficult and yet so fundamental for many 
branches of science that each concrete result is useful in shedding light on 
the nature of general solution. 

In spite of some bogus claims [55] , refuted in [35] , this result had stood for more 
then three decades as the only known solution of the iV-representability prob- 
lem beyond two electrons A 2 H r and two holes A r ~ 2 H r . For the latter systems 
the problem is easy and the constraints amounts to double degeneracy of the 
spectrum, starting from the head \2i-1 — A 2 i for two electrons and from the tail 
A r -2i = A r _2i-i for two holes [5], where we set Ai = for i > r, and Ai = 1 for 
I < 1. 

Here we solve this longstanding problem. The content of the paper is as 
follows. 

In Section[5Jwe recast the Berenstein-Sjamaar theorem [T] Thm 3.2.1] into a 
usable form (Theorem [T]) . This provides a theoretical basis for our study. 

We start Section[3jby a variation of the above problem, called v -representabi- 
lity^ that takes into account both spin and orbital occupation numbers. Mathe- 
matically this amounts to replacing the exterior power f\ N TL by a representation 
TL U defined by Young diagram v of order N. Theorem [2] gives a formal solution 
of the I'-representability problem. We derive from it the majorization inequality 
A ^ is, that plays the role of the Pauli principle. This inequality is necessary and 
sufficient for A to be occupation numbers of an unspecified mixed state (Theorem 
[3j). Theorem [4] deals with a class of systems where the majorization inequality 
alone provides a criterion for pure i/-representability. This includes the so-called 
closed shell, meaning a system of electrons of total spin zero. The corresponding 
Young diagram v consists of two columns of equal length. For this system all 
constraints on the occupation numbers are given by the Pauli type inequality 
A < 2. In the next Theorem [5] we calculate the topological coefficients c v w (a) that 
governed the constraints on the occupation numbers in Theorem [2] This gives 
it the full strength we need in the next section. 

Section [¥] starts with analysis of pure i^-representability for a toy example of 
two- row diagrams, that allows us to illustrate the basic technique (Theorem [6]) . 

2 Recently M.B. Ruskai published the proof 33_ derived from known constraints on the spec- 
tra of Hermitian matrices A, B, and C = A + B. Conceptually the 7V-representability problem 
is close to the Hermitian spectral problem 15 16], but a direct connection between them, be- 
yond sporadic coincidences, is unlikely. An independent R.L. Kingsley's solution apparently 
has never been published. 
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These are exceptional systems where the constraints on the occupation numbers 
are given by a finite set of inequalities independent of the rank. Then we return 
to the original ./V-rcprcscntability problem, that appears to be the most difficult 
one. For example, in contrast to Theorem [6l no finite system of inequalities can 
describe 7V-rcprcscntability for a fixed N > 1 and arbitrary big rank (Corollary 
|3]to Proposition [5} • This forces us to restrict either the rank, as we do in the 
last section, or the type of the inequalities. Here we focus on the inequalities 
with 0/1 coefficients. It turns out that under some natural conditions such an 
inequality should be either of the form 

Ai x +K 2 +--- + A ijv _ 1 < N-2, (2) 

with X)fc(*fc — k) = T — N + 1, or of the form 

A il +A ia + --- + A ip <JV-l, (3) 

with p> N and X)fc(*fc — k) = We call them Grassmann inequalities of the 
first and second kind respectively. A surprising result is that these inequalities 
actually hold true with very few exceptions (Theorems [7] and [SJ . 
In the simplest case N = 3 we get from ([5} inequalities 

A fc+ i + A r _fe<l, 0<fc<(r-l)/2 

that hold for any even rank r > 6. This constraint prohibits more than one 
electron to occupy two symmetric orbitals and supersedes the original Pauli 
principle. For r = 6, due to the normalization J^. \ = 3, the inequalities degen- 
erate into Borland-Dennis equalities ([1]) . For odd rank the first inequality k = 
should be either skipped or replaced by weaker one Ai + A r < 1 + — ^j. 

We treat Grassmann inequalities of the second kind ([3]) only for lowest levels 
p = N, N + 1 . For N — 3 and p — N + 1 they amount to four inequalities 

A 2 + A 3 + A 4 + A 5 < 2, Ai + A 3 + A 4 + A 6 < 2, 

Ai + A 2 + A 5 + A 6 < 2, Ai + A 2 + A 4 + A 7 < 2, w 

that hold for arbitrary rank r and give all the constraints for r < 7. For r = 6 
they turn into Borland-Dennis conditions (fTJ). 

In the next Section [5] we briefly discuss a connection of the t'-representability 
with representation theory, that provides information complementary to Theo- 
rem [2] A combination of the two approaches leads to an algorithm for solution 
of the problem for any fixed rank. The algorithm, along with other tools, has 
been used in calculations reported in the last Section [6l Eventually this led to a 
complete solution of the ./V-representability problem for rank r < 10. However, 
we provide a rigorous justification only for r < 8. We also give an example of 
constraints on the spin and orbital occupation numbers for a system of three 
electrons of total spin 1/2. 

The first sections may be mathematically more demanding then the rest of the 
paper. We recommend books [ZUHIE] as a general reference on Schubert calculus, 
Lie algebra, and representation theory. 

The theoretical results of the paper belong to the second author. They were 
often inspired by calculations, that at this stage couldn't be accomplished by a 
computer without intelligent human assistance and insight. 
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2. A digest of the Berenstein-Sjamaar paper 

Let M be a compact connected Lie group with the Lie algebra m and its dual 
coadjoint representation m*. For coadjoint orbit O C m* of group M and a 
Cartan subalgebra t C m consider the composition A : O <—>■ m* — > t* known as 
the moment map. By Kostant's theorem its image is a convex polytope spanned 
by the IF-orbit of some weight // G t* which can be taken from a fixed positive 
Weyl chamber t+. Here W — N(t)/Z(t) is the Weyl group of M. This gives a 
parameterization of the coadjoint orbits by the dominant weights fi G t?_. 

Example 1. In this paper we will mostly deal with the unitary group U(n) whose 
Lie algebra u(n) consists of all Hermitiaifl n x n matrices. Let us identify u(n) 
with its dual via the invariant trace form (A, B) = Tr(AB). Then the (co)adjoint 
orbit consists of all Hermitian matrices A of spectrum /i : fi± > \i% > ■ ■ • > fi n 
and the moment map A : O ^ — ^ t is given by orthogonal projection into the 
Cartan subalgebra of diagonal matrices t. Kostant's theorem in this case amounts 
to Horn's observation that the diagonal entries of Hermitian matrices of spectrum 
[i form a convex polytope with vertices u>/i obtained from /i by permutations of 
the coordinates This is equivalent to the majorization inequalities 

di < /ia 

^1+^2+^3 < Ml + ^2 + (5) 



di + d 2 H h d n = /ii + A*2 H h Mn 

for the diagonal entries d : d\ > d% > • • • > d n of matrix ^4. We will use for them 
a shortcut d ^ ji. 

Consider now an immersion / : L — * M of another compact Lie group L and 
the induced morphisms /*:['—> m and /* : m* — > C* of the Lie algebras and their 
duals. In the paper pQ Berenstein and Sjamaar found a decomposition of the pro- 
jection C 1* of an M-orbit O p C m* into i-orbits O x C Here we 
paraphrase their main result in the form suitable for the intended applications. 

Fix a Cartan subalgebras II iu and for every test spectrum a G tz, consider 
the inclusion of the adjoint orbits of groups L and M 

Va :O a ^O fAa) (6) 

through a and /* (a) respectively. Topologically the orbits are (generalized) flag 
varieties. They carry a hidden complex structure coming from the representation 

O a = L/Z L (a)=L c /P a (7) 

where P a C L is a parabolic subgroup of the complexified group L c whose Lie 
algebra p a is spanned by tt, and the root vectors X a such that (a, a) > 0. One 
can say this in another way 

P a = {g G L c | lim e* Q 5e-* a exists} 



3 Hereafter we treat u(n) as the algebra of Hermitian, rather than skew-Hermitian, operators 
at the expense of a modified Lie bracket [X, Y] = i(XY — YX). 
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which makes it clear that / : P a — > Pf t ( a ) . 

We will use the parabolic subgroups to construct canonical bases in cohomolo- 
gies H*(O a ) and H*(Of t ^). Let T L C B C P a be a Borel subgroup containing 
a maximal torus Tj, with Lie algebra t/,. The flag variety O a — L c /P a splits 
into disjoint union of Schubert cells BvP a /P a , parameterized by the left cosets 
v £ Wl /Wz L ( a ) or in practice by representatives of minimal length I = £(v) 
in these cosets. We actually prefer to deal with shifted cells v~ 1 BvP a /P a = 
B v P a /P a depending on the Borel subgroups B v D Tl modulo conjugation by 
the Weyl group of the centralizer W (Z z,(a)) . The closure of B v P a /P a is known 
as the Schubert variety, and its cohomology class a v £ H 2l ^ v \O a ) is called the 
Schubert cocycle. These cocycles form the canonical basis of the cohomology ring 
H*(O a ). 

Inclusion ([6]) induces a morphism of the cohomologies 

<p* a :H*(0 Ma) )^H*(O a ), (8) 
given in the canonical bases by the coefficients c^, (a) of the decomposition 

<P*a °w i-> ^c v w (a)a v . (9) 

V 

They play a crucial role in the next theorem. We extend them by zeros if either 
v £ Wl or w £ Wm is not the minimal representative of a coset in Wl/Wz^u) 
or W M /W ZM ( fit ( a)) respectively. 

Theorem 1. In the above notations the inclusion 0\ C is equivalent to 

the following system of linear inequalities 

(A, va) < (fi,wf*(a)) (a,v,w) 

for all a e i L , v <G Wl,w £ Wm such that c v w (a) ^ 0. 

Proof. This is not the way how Berenstein and Sjamaar stated their result. 
Instead, for some generic oo 6 U they fix positive Weyl chambers 9 oq 
and t M 3 /*(a ) and use them to define Schubert cocycles a v £ H*(O a ) and 
o w £ H*{Of_^ a )) for all other a £ tj. Hence their Schubert cocycles a w are 
canonical in the above sense iff /*(a) and /*(ao) are in the same Weyl chamber. 
The set of such a £ tj form a convex polyhedral cone called the principle 
cubicle. It is determined by a$, and different choices of ao produce a polyhedral 
decomposition of the positive Weyl chamber into cubicles. 

For every cubicle Berenstein and Sjamaar gave a system of linear constraints 
on the dominant weights A, /i, so that all together they provide a criterion for 
the inclusion 0\ C For the principal cubicle the constraints are most 

simple and look as follows [U Thm 3.2.1] 

rtef^V-C), for 4(a )^0, (10) 

where C is a cone spanned by the positive roots in t* M . Note that f*(C) is the 
cone dual to the principal cubicle and therefore the above condition can be recast 
into the inequalities 



(w -1 A,o) < <r(urV),a) (A,«o) < (n,u>Ma)), (11) 
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that hold for all a from the principle cubicle provided that c^(ao) 0- The 
coefficients c^(a) are actually constant inside the cubicle, and therefore the last 
condition can be changed to c^(a) ^= 0. Thus we arrived at the inequalities 
( jq, v , for the principle cubicle. Other inequalities ( ja, v, mj ) follow by choosing 
another cubicle as the principle one. They are equivalent to the remaining more 
complicated inequalities in [IJ Thm 3.2.1], but look differently since Berenstein 
and Sjamaar use other non-canonical Schubert cocycles. □ 

Example 2. Quantum marginal problem \17\ . Let's illustrate the above theorem 
with immersion of unitary groups 

/ : V(H A ) x \J{H B ) -» U(«ab), 9a^9b^ 9a ® 9b, 

where TCab — 7~(-A 8%. As we have seen in Example [T] the coadjoint orbit of 
U(Wab) consists of the isospectral Hermitian operators pab ■ 7~Lab understood 
here as mixed states. The projection 

}*{pab) — Pa ® 1 + 1 ® Pb 

amounts to reduced operators pa ■ Ti-A and ps '■ Ti-B implicitly defined by the 
equations 

^ Ha {paXa) = Ti nAB (p AB X A ), Tt Hb {paX b ) = Ty Hab (pabX b ) (12) 

for all Hermitian operators Xa ■ Ti-A and Xb '■ T~Cb- This means that pa, Pb are 
just the visible states of the subsystems Ha, 1~Lb- In this settings Theorem Q] tells 
that all constraints on the decreasing spectra \ AB = Spec(p AB ), X A = Spec(p A ), 
and X B = Spec(p B ) are given by the inequalities 

E a * A «« + E b ^u) ^ E(° + h )i x tf k) , (is) 

i j k 

for all test spectra a : a\ > 02 > ■ ■ ■ > a n , b : b\ > b2 > ■ ■ ■ > b m from the 
Cartan subalgebras tA, ts and permutations u, v, w such that c™(a, b) 7^ 0. Here 
(a + by denotes the sequence a, + &j arranged in decreasing order. The order 
determines the canonical Weyl chamber containing /*(<z, b). The pairs (a, 6) with 
fixed order of terms at + bj in (a + b)^ form a cubicle. 

The adjoint orbit O a C u(H.a) is a classical flag variety understood as the 
set of Hermitian operators Xa ' Ti-A of spectrum a = Spec X a- Denote it by 
FaiH a)- Then the morphism ((6J is given by the equation 

Vab : Ta{U A ) x T b (H B ) -» F a+b {H A B), {X A , X b ) ^ X a ®1 + I ® X b (14) 

and the coefficients c™(a, b) are determined by the induced morphism of the 
cohomologies 

<p* ab : H*(T a+b (H AB )) -> H*(T a (H A )) <8 H*(T b (H B )) 

<r w ^^2cl v (a,b) ■ a u (g> a v . (15) 

One can find the details of their calculation in [T7]. Note that c™(a, 6) = 1 for 
identical permutations u, u, io. Hence we get for free the following basic inequality 

E aiXf + E < J> + 6)iA^ s (16) 

i j fc 

valid for all test spectra a, 6. 
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3. One point ^-representability 

In this section we apply the above results to the morphism / : U(7Y) — > U(7Y !/ ) 
given by an irreducible representation Tt v of group \j(Ti) with a Young diagram 
v of order N = \v\. For a column diagram we return to the A^-fermion system 
A N H, while a row diagram corresponds to the iV-boson space S N TL. However, 
the main reason to consider the general para- statistical representations TL V is 
not a uniform treatment of fermions and bosons, but taking into account spin. 
Observe that the state space of a single particle with spin splits into the tensor 
product H — H r ®H. S of the orbital H r and the spin Tt s degrees of freedom. 
The total iV-fermion space decomposes into spin-orbital components as follows 

m 

/\ N {Hr ®H S )= Ys H r® K* > ( 17 ) 

\v\=N 

where v stands for the transpose diagram. In many physical systems, like elec- 
trons in an atom or a molecule, the total spin is a well defined quantity that 
singles out a specific component of this decomposition. Theorem [1] applied to 
the component gives all constraints on the possible spin and orbital occupation 
numbers, see the details in n° l3.1.1l below. 



3.1. Physical interpretation. Let's now relate Theorem 1 to the iV-representabi- 
lity problem and its ramifications indicated above. We'll refer to the latter as 
the v-representability problem. 

It is instructive to think about X £ u(7Y) as an observable and treat p S u(7i)* 
as a mixed state with the duality pairing given by the expectation value of X in 
state p 

(X,p)=Tr n Xp (18) 

(forget for a while about the positivity p > and normalization Trp = 1). 

We want to elucidate the physical meaning of the projection /* : u(W)* — * 
u(H)* uniquely determined by the equation 

(f*(X), P »} = {X,f*(p v )}, Xeu(H), p"eu(?0*. 

In the above setting (fl8|) it reads as follows 

Tr- H "(Xp u ) = r KH(Xr(p v )), VXeu(H). (19) 

A good point to start with is Schur's duality between irreducible representations 
of the unitary Xj(7i.) and the symmetric Sn groups 

H m = K v ®S V . (20) 

\v\=N 

The latter group acts on 7-t® N by permutations of the tensor factors, and its 
irreducible representations S u show up in the right hand side. One can treat 
7Y®JV ag a s t a t e S p aC e of A^-particles, and for identical particles all physical 
quantities should commute with Sn. Looking into the right hand side of (|20(1 we 
see that such quantities are linear combinations of operators p v ® 1 acting in the 
component W ® S v and equal to zero elsewhere. In the case of a genuine mixed 
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state p v ', i.e. a nonnegative operator of trace 1, one can treat (p v ® 1)/ d\mS v 
as a mixed state of N identical particles obeying some para-statistics of type v. 
Let pi : Ti. be its i-th reduced state. Since p v £g> 1 commutes with Sn, the reduced 
state p — pi is actually independent of i. However, occasionally we retain the 
index i just to indicate the tensor component where it operates. 

Proposition 1. In the above notations 

f*(p u )=Np. (21) 
Proof. We have to check that ([21]) fits the equation (fT9|) : 

i 

where Xi is a copy of X acting in the i-th component of H® N , so that 

p v <g> 1 

Tr w ®« A, — — = Tr w Ajpj 

by definition (fT!Z|) of reduced state. □ 

A general v-representability problem concerns with the relationship between 
the spectrum p, of a mixed state p" and spectrum A of its particle density matrix 
N p. The latter spectrum is known as the occupation numbers |3 of the system in 
state p v . Formally the constraints on the spectra are given by Theorem [TJ 

Remark 1. The above construction allows for a given mixed state p v to define 
the higher order reduced matrices. Their characterization would have almost 
unlimited applications. Indeed, behavior of most systems of physical interest is 
governed by two-particle interaction. As a result, the energy of a state becomes 
a linear functional of its two-point reduced matrix. To minimize the energy and 
to find the correlation matrix of the ground state one has to elucidate all the 
constraints that a two-point reduced matrix should satisfy. This problem and the 
whole program are known as Coulson challenge^ [6]. In the form just described 
it may be unfeasible even for quantum computers [23]. For other approaches 
and the current state of art see [26] . This problem is far beyond the scope of our 
paper. Nevertheless, the characterization of one point reduced matrices given 
below imposes also new constraints on the higher reduced states. 



3.1.1. Constraints on spin and orbital occupation numbers. Let's return to a 
system of N fermions, this time of smallest possible spin s = 1/2, dim7i s = 2. 
In this case spin-orbital decomposition (fT7| involves only terms 

U v r <g> W* (22) 

with at most two-column diagram v. The sizes of the columns a > (3 are deter- 
mined by equations 

a + fJ = N, a -(3 = 2 J, (23) 

4 More precisely, the occupation numbers of natural orbitals. The latter are defined as eigen- 
vectors of the particle density matrix. 

5 And also as two-particle ./V-representability or, following D. Herschbach, as a holy grail of 
theoretical chemistry. 



The Pauli principle revisited 



9 



where J is the total spin of the system, so that Ti. v s — 7ij is just the spin J 
representation of the group S\J(H S ) = SU(2). 

Consider now a pure iV-fermion state of total spin J 

where the diagram v is determined by equations (|23[) . Let p v and p J be its 
reduced states in the orbital and spin components respectively. The basic fact is 
that the reduced states are isospectral Spec p v — Spec p J . Hence Spec p v can be 
identified with the spin occupation numbers. On the other hand Theorem [TJ in 
view of Proposition [TJ relates Specp^ with the orbital occupation numbers given 
by the spectrum of the particle density matrix Np. In this way one can produce 
all constraints on allowed spin and orbital occupation numbers, provided that 
a solution of the i^-representability problem is known for two-column diagrams. 
We address this issue in sections [3T21 and [3731 See also Corollary [JJ in section l3~2l 

3.2. Formal solution of the v-representability problem. Henceforth we treat the 
lower index r as the rank of the Hilbert space Tt r . Recall that the character of 
the representation TL", i.e. the trace of a diagonal operator 

Z = AlB,g{zx,Z2, ...,z r ) 6U(H r ), (24) 

in some orthonormal basis e of Tl r , is given by Schur's function S v (zi, z%, . . . , z r ). 
It has a purely combinatorial description in terms of the so called semistandard 
tableaux T of shape v. The latter are obtained from the diagram v by filling it 
with numbers 1, 2, . . . , r strictly increasing in columns and weakly in rows. Then 
the Schur function can be written as a sum of monomials z T = JligT z i 

sM = Y. zT 

T 

corresponding to all semistandard tableaux T of shape v. The monomials are 
actually the weights of representation Ji", meaning that 

z ■ ex = z T er (25) 

for some basis ex of l-i v r parameterized by the semistandard tableaux. Denote 
by t C u(H r ) and C u(H^) the Cartan subalgebras of real diagonal operators 
in the bases e and ex respectively, so that the differential of the above group 
action z : ct i— * z T ex gives the morphism 

/* : t -> lu, /»(o) : e T i-> a T e T , (26) 

where ax '■= X)ieT a «- ^ s m Example [2] we treat the orbits O a and C/»( a ) as 
flag varieties T a {H r ) and T^iTi-r) consisting of Hermitian operators of spectra 
a : <2i > ct2 > • • • > a r and a v respectively. Here a v consists of the quantities olt 
arranged in the non-increasing order 

a v := {ar \ T — semistandard tableau of shape v}K (27) 

Finally, we need the morphism 



(28) 
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(29) 



(30) 



V 



Theorem 2. In the above notations all constraints on the occupation numbers 
A of the system l~L v r in a state p v of spectrum [i are given by the inequalities 



for all test spectra a and permutations v,w such that c^(a) 0. 

Proof. In view of Proposition[Tl this is what Theorem[T]tells. One has to remem- 
ber that the left action of a permutation on "places" is inverse to its right action 
on indices. That is why the permutations v and w, acting on a and f*(a) = a v 
in Theorem [TJ move to the indices of A and /i in the inequality ([3Tjl . □ 

The coefficient c^ (a) depends only on the order in which quantities or appear 
in the spectrum a". The order changes when the test spectrum a crosses a 
hyperplane 



The hyperplanes cut the set of all test spectra into a finite number of polyhedral 
cones called cubicles. For each cubicle one has to check the inequality (f3Tj) only 
for its extremal edges. As a result, the z^-representability amounts to a finite 
system of linear inequalities. 

Remark 2. Let's emphasize once again the difference between Bcrenstein-Sjamaar 
theorem p] Thm 3.2.1] and its version used in this paper. In the settings of The- 
orem [3] it manifests itself in the way how the quantities ax are ordered in the 
spectrum a", or what is the same which parabolic subgroup is used for definition 
of Schubert cocycles. Berenstein and Sjamaar choose a specific order of tableaux 
T, while we rely on the natural order of the quantities ax = Sigt The letter 
choice allows to treat the inequalities uniformly, and to avoid a rather cumber- 
some transformation every time the test spectrum passes from one cubicle to 
another. 

Recall from n° 13.1.11 that the theorem also describes a relationship between 
the spin and orbital occupation numbers. We keep for them the above notations 
fi and A respectively. 

Corollary 1. All constraints on spin and orbital occupation numbers of N- 
electron system in a pure state of total spin J are given by the inequalities (|31|) , 
applied to two column diagram v determined by equations (|23p . and bounded to 
mixed states p v of rank not exceeding dimensionality 2 J + 1 of the spin space. 

a 

We postpone the calculation of the coefficients c^(a) to section 13731 and focus 
instead on some general results that can be deduced from the theorem as it 
stands. 




(31) 



k 
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Being a ring homomorphism, ip* maps unit into unit 
1 for identical permutations v,w. Hence the following 

i k 

holds for all test spectra a. Let's look at it more closely for a pure state p v = 
\ip)(tp\ in which case the right hand side is maximal and the inequality takes the 
form 

cijXi < a" = maxT^ a, = am, (32) 

i i£T i 

where v\ > V2 > ■ ■ ■ > are rows of v. The maximum in the right hand side is 
attained for the tableau T of shape v whose i-row is filled by i. 

The normalization Xi = N — yV Vj allows to shift the test spectra into 
the positive domain a\ > a 2 > • • • > 0, so that they became nonnegative linear 
combinations of the fundamental weights 

u k = (1,1,. ..,1,0,0,. ..,0). (33) 
fe 

Hence it is enough to check (|32[) for a = lju, that gives the majorization inequality 
X < v : cf. Example [T] Thus we arrived at the first claim of the following result 
that characterizes occupation numbers of system 7i v in an unspecified mixed 
state. 

Theorem 3. The occupation numbers of the system TL U in an arbitrary mixed 
state satisfy the majorization inequality 

X<v, (34) 

and any such X can be realized as the occupation numbers of some mixed state. 

Proof. The second claim follows from two observations: 

1. The occupation numbers of a coherent state ip € TL V , that is a highest vector 
of the representation, are equal to v. 

2. The set of allowed occupation numbers, written in any order, form a convex 
set. 

Indeed, the polytope given by the majorization inequality (|34j) is just a convex 
hull of vectors obtained from v by permutations of coordinates, cf. Example [TJ 
Hence by 1 and 2 it consists of legitimate occupation numbers. 

Proof of 1. Consider a decomposition of the complexified Lie algebra 

u(H) ® C = Ql(H) = n_ + fj + n +! 

into a diagonal Cartan subalgebra f) = t® C accompanied with lower- and upper- 
triangular nilpotent subalgebras n^. By definition n+ annihilates the highest 
vector ip £ U v of weight v. Hence (ip\X ± \ij) = {X^if>\ip) = for all X ± £ n±. 
Then by equation IT9|) 



3.2.1. Basic inequalities. 
<p* a {l) = 1, that is t*(o) = 
basic inequality 



(vi*^) = Tr^p^x^i) - Tv n (x ± r(W(m = o, yx ± g n±. 
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This means that p = f*{\ip){ij)\) is a diagonal matrix. On the other hand tip = 
(t, v)ip for t G t, hence as above 

(t,u) = ^\t\if,) = TY W »(*MM) =Tr n (trm(m =Tr H (tp) = (t,p), 
that is Spec p = v. 

Proof of 2. Let p\, p^ be mixed states, with the particle densities pi, p2, and 
the occupation numbers Ai, A2. Apply to pi, p\ a unitary rotation pi 1— > UpiU*, 
p\ 1— ► U p"U* that transforms orthonormal eigenvectors of pi into that of pi in 
a prescribed order. The resulting new operators p±, pi commute and have the 
original spectra Ai, A2. Then the particle density matrix p = pipi + P2P2 of the 
convex combination p v = pip\ + P2P2 nas spectrum A = piXi + P2^2- □ 

For a column diagram v the majorization inequality A -< v amounts to the 
Pauli exclusion principle Xi < 1. In general, we refer to it as the Pauli constraint. 
Note that the above proof shows that equality in (|34f is attained for the coherent 
states only. The second part of Thcorem[3]extends Coleman's result !5 for h N TL. 

Recall, that the theorem solves the ^-representability problem for unspecified 
mixed states. We will see later that for pure states the answer in general is much 
more complicate. Nevertheless, there are surprisingly many systems for which 
the majorization inequality along is sufficient for pure ^-representability. We 
address them in the next item. 

3.2.2. Pure moment polytope. One of the most striking features of Theorem [2] is 
the linearity of the constraints ([31]) . As a result, the allowed spectra (A, p) form a 
convex polytope, called (noncommutative) moment polytope. The convexity still 
holds for any fixed p = Spec p v , and in particular for the occupation numbers 
A of all pure states. We refer to the latter case as the pure moment polytope. It 
sits inside the positive Weyl chamber, and its multiple kaleidoscopic reflections 
in the walls of the chamber generally form a nonconvex rosette, consisting of all 
legitimate occupation numbers written in an arbitrary order. It can be convex 
only if all constraints on the occupation numbers are given by the majorization 
inequality \ < v alone. Here we describe a class of representations H u with this 
property. 

This happens, for example, for a system of N > 2 bosons. In this case v is 
a row diagram and the majorization inequality imposes no constraints on A. By 
Theorem [3] this means that every nonnegative spectrum A of trace N represents 
occupation numbers of some mixed state. However for bosons one can easily find 
a pure state that does the job: 

* = E y/** e ? N e sN ^ 

i 

where is an orthonormal basis of Tt. This makes the bosonic Af-representability 
problem meaningless. 

A more interesting physical example constitutes the so-called closed shell, 
meaning a system of electrons of total spin zero. The corresponding diagram v 
consists of two columns of equal length. We will see shortly that in this case the 
Pauli constraint A < 2 shapes the pure moment polytope. 

Observe that it is enough to construct pure states whose occupation numbers 
are generators of the cone cut out of the Weyl chamber by the majorization 
inequality \ <v. Then the convexity does the rest. 
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Recall, that in the proof of Theorem [3] we have already identified v with 
the occupation numbers of a coherent state. Due to the majorization inequality 
A ^ v, the entropy of its reduced state is minimal possible. By that reason 
coherent states are generally considered as closest to classical ones [30]. At the 
other extreme one finds the so called completely entangled states ip € H." whose 
particle density matrix p = f*{\if})(ip\) is scalar and the reduced entropy is 
maximal [19]. By definition (OH we have Tr H (Xp) = Tr n ,{X\t/j)(ip\) = (ip\X\tp), 
so that the completely entangled states can be described by equation 

(ip\X\ip)=Q, Vlesu(H). (35) 

Let's call a system "H" exceptional if the SU(7i r ^representation H." is equivalent 
to one of the following: H r , its dual H*, and, for odd rank r, A 2 H r , /\ 2 H*. The 
Young diagram v of an exceptional system can be obtained from r x m rectangle 
by adding an extra column of length 1, r — 1,2, r — 2 respectively. 

One readily realizes that the exceptional systems contain no completely en- 
tangled states, say because reduced matrix of tp € f\ 2 TL r has an even rank. 

Proposition 2. In every non- exceptional system TL V there exists a completely 
entangled state. 

Proof. The result is actually well known, but in a different context. The entan- 
glement equation (|35[) is nothing but the stationarity condition for the length 
of vector (ip\ip) with respect to action of the complexified group SL(H). It is 
known |34j that every stationary point is actually a minimum, and an SL(7i)- 
orbit contains a minimal vector if and only if the orbit is closed. As a result, we 
end up with the problem of existence of a nonzero closed orbit, or, what is the 
same, the existence of a nonconstant polynomial invariant. The proposition just 
reproduces a known answer to the latter question 34 . □ 

By admitting other simple Lie groups we find only two more exceptional 
representations: the standard representation of the symplectic group Sp(n) and 
a halfspinor representation of Spin(10). 

Now we can solve the pure i^-representability problem for a wide class of 
systems, including the above mentioned closed shell. 

Theorem 4. Suppose that all columns of Young diagram v are multiple, meaning 
that every number in the sequence of columns lengths v\ > v\ > v\ > • • • appears 
at least twice. Then all constraints on the occupation numbers of the system Ti v 
in a pure state are given by the majorization inequality A -< v along. 

Proof. We'll proceed by induction on the height of the diagram v. The trivial- 
ity of the bosonic TV-representability problem provides a starting point for the 
induction. 

Let now A be a vertex of the polytope cut out of the positive Weyl chamber by 
the majorization inequality \ <v. Take notice that the latter includes equation 
Tr A = Trv. Then the following alternative holds: 

1 . Either all nonzero components of A are equal, 

2. Or one can split A and v into two parts A = A'|A", v = v'\v" containing the 
first p components and the remaining ones, both satisfying the inequalities 
A' r< u', A" r< v" . 
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Indeed, the second claim just tells that the p-th majorization inequality in |5| 
turns into equation. On the other hand, if all the majorization inequalities are 
strict, and A contains different nonzero entries, then one can linearly vary these 
entries preserving the non-increasing order of A and the majorization A ^ v. As 
result we get a line segment in the polytope containing A, which is impossible 
for a vertex. 

We've to prove that every vertex A represents occupation numbers of some 
pure state. Consider the above two cases separately. 

Case 1. Let A contains r equal nonzero entries and H r C H be a subspace of 
dimension r. The conditions of the theorem ensure that the system H v r is non- 
exceptional, hence by Proposition [2] it contains a state ^ € H v r with occupation 
numbers equal to nonzero part of A. In bigger system H u D H v r its occupation 
numbers will be extended by zeros. 

Case 2. Let the system has rank r = p+q. Choose a decomposition H r = H p ®H q 
and consider a restriction of the representation H v r onto subgroup U(H P ) x XJ(H q ) 

where c" are the omnipresent Littlewood-Richardson coefficients. Observe that 
(fyiyii = 1, and therefore Hp <g> H q C H v r . By induction hypothesis there exist 
states ip' S Hp and if>" € H v q with occupation numbers A', A" and particle 
densities p' , p" respectively. Then decomposable state i/j = ip' <g> ip" has particle 
density p' © p" , and its occupation numbers are equal to A = A'|A". □ 

Let's extract for a reference a useful corollary from the last part of the proof. 

Corollary 2. Suppose that the Littlewood-Richardson coefficient is nonzero. 
Then merging of the occupation numbers A', A" of the systems Hip 1 , H q form 
legitimate occupation numbers of the system Hp +q . □ 

Remark 3. The restriction on the column's multiplicities of the diagram v is 
needed only to ensure that the components of any splitting v = v'\v"\v"'\ . . . are 
non-exceptional. The latter condition holds for any two-row diagram [a,/3],/3 7^ 
1 for dim7i > 3. This gives examples of systems beyond Theorem [4] say for 
v = [3, 2], whose pure moment polytope is given by the majorization inequality 
along. More such diagrams can be produced as follows: take v as in Theorem |4] 
and remove one cell from its last row. This works when the last row contains at 
least three cells and rank of the system is bigger than the height of v. A complete 
classification of all such systems is still missing. 

3.2.3. Dadok-Kac construction. In the last two theorems we encounter the prob- 
lem of construction a pure state with given occupation numbers. The problem 
lies at the very heart of the z/-representability and one shouldn't expect an easy 
solution. Nevertheless, there is a combinatorial construction that produces a 
state with diagonal density matrix, whose spectrum can be easily controlled. 
It has been used first by Borland and Dennis [3] to forecast the structure of 
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the moment polytope for small fermionic systems. Later on Miiller [27| formal- 
ized and advanced their approach to the limit. It fits into a general Dadok-Kac 
construction [lOj that works for any representation. 

Below we follow the notations introduced at the beginning of n° 13.21 Let 
x = diag(a;i, X2, ■ ■ ■ , x r ) be a typical element from Cartan subalgebra t C u(7i r ). 
For a given semi-standard tableau T call the linear form lot ■ % %t = X^eT x i 
the weight of the basic vector ct G H-r- We also need nonzero weights of the 
adjoint representation Oiq^j • X ' ' Xi Xj , i 7^ j called roots. Let's turn the set 
of semi-standard tableaux of shape v into a graph by connecting T and T' each 
time lot — is a root, i.e. the contents of T and T", considered as multi-sets, 
differ by exactly one element. 

Proposition 3. Let T be a set of semi- standard tableaux of shape v containing 
no connected pairs. Then every state ip = Xtst °t&t € W with support T has 
a diagonal particle density matrix with entries 



where every tableau T is counted as many times as the index i appears in it. 

Proof. The proof refines the arguments used in claim 1 of Theorem[3l from which 
we borrow the notations. As in the above theorem we have to prove (ip\X\ip) = 
for every X 6 n+ + n_. It is enough to consider root vectors X a that form a 
basis of n + + n_ . Then 



Since X a eT has weight a + lot, it is orthogonal to ey, except lot 1 = lot + ol. 
The latter is impossible for T, T' S T, and therefore the reduced state of ip is 
diagonal. A straightforward calculation gives the diagonal entries ([37)) . □ 
We'll have a chance to use this construction in n° 14.11 

Take notice that for a fixed support T the set of unordered spectra (|37|) 
form a convex polytope. It is not known when this approach exhausts the whole 
moment polytope. The smallest fermionic system where it fails is /\ 3 Hs, see n°[6l 

3.3. Calculation of the coefficients c^(a). To move further and to give Theorem 
Q]thc full strength one has to calculate the coefficients c^(a). Berenstein and 
Sjamaar left this problem mostly untouched. However, in the i^-representability 
settings, highlighted in Theorem [21 this can be done pretty explicitly. 

3.3.1. Canonical generators. To proceed we first need an alternative description 
of the cohomology of flag variety T a (Ji-r) [Zj- Recall that the latter understood 
here as the set of Hcrmitian operators in TL r of given spectrum a. To avoid 
technicalities, we assume the spectrum to be simple a\ > 02 > ■ ■ ■ > a r . Let 
Si be the eigenbundle on JF a (?i r ) whose fiber at X € T a (7i r ) is the eigenspace 
of operator X with eigenvalue Oj. Their Chern classes Xi — Ci(£i) generate the 
cohomology ring H*{T a {l-L r )) and we refer to them as the canonical generators. 
The elementary symmetric functions o~i (x) of the canonical generators are the 




(37) 



{iP\X a \i ) )= CT>c T (eT'\X a \e T ). 



T,T'eT 
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characteristic classes of the trivial bundle TL r and thus vanish. This identifies 
the cohomology with the ring of coinvariants 

H*(J r a (H r )) = Z[xi,x 2 , . . . ,x r ]/ (at, er 2 , . . . ,oy). (38) 

This approach to the cohomology is more functorial and by that reason leads to 
an easy calculation of the morphism (|29p 

<p* a : H*(F a »(H v )) H*(T a (H)). 

Recall that the spectrum a v consists of the quantities or = J2i£T ai arranged in 
decreasing order, where T runs over all semi-standard tableaux of shape v. We 
define xt = SieT x i m a similar way. 

Proposition 4. Let Xi and x k be the canonical generators of H*(J- a (Tl)) and 
H* (J^a* (?i v )) respectively. Then 

tPa( x k) = x t, when a k — ax- (39) 
In other words, ip*(x^.) is obtained from a k by the substitution ai <— > Xi. 

Proof. The eigenbundle £i is equivariant with respect to the adjoint action X i— > 
uXu* of the unitary group \](Ti.). Therefore it is uniquely determined by the 
linear representation of the centralizcr D = Z(X) in a fixed fiber £i(X) or 
by its character Si : D — > S 1 = {z e C* | \z\ — 1}. In the eigenbasis e of 
the operator X the centralizer becomes a diagonal torus with typical element 
z = diag(zi, Z2, ■ ■ ■ , z r ) and the character s\ : z i— > Zj. 

Let now X" = <p a (X), D v = Z(X V ), and e T be the weight basis of H' y , 
introduced in section 13 . 21 parameterized by semi-standard tableaux T of shape 
v and arranged in the order of eigenvalues a". Then the character of the pull 
back <p~ 1 {£ k ) is just the weight JlieT e « 01 the k-th vec t° r 6t, where the tableau 
T is determined from the equation a k = ar, cf. (|25p. Thus < J 9~ 1 (£^) = (££>j gT £j 
and we finally get 

<p*M) = = cifaHw) = c!(0^) = 53^ = ^. □ 

Remark 4- Formula (|39[) may look ambiguous for a degenerate spectrum a, while 
in fact it is perfectly self-consistent. Indeed, consider a small perturbation a, 
resolving multiple components of a, and the natural projection 

that maps X = J^. | e^} {e^ | into X = J^. aj|ei)(ej|, where is an orthonormal 
eigenbasis of X. It is known [2] that 7r induces isomorphism 

7T* : H*{T a (H)) ~ H*(F~ a {H)) w ( D \ (40) 

where on the right hand side stands algebra of invariants with respect to per- 
mutations of the canonical generators ii with the same unperturbed eigen- 
value ai — a. Such permutations form Weyl group W(D) of the maximal torus 
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D = Z(X) in D = Z(X). For example, characteristic classes of the eigenbun- 
dlc £ a with multiple eigenvalue a = correspond to elementary symmetric 
functions of the respective variables x i . 

Equation (|39|) . as it stands, depends on a specific ordering of the unresolved 
spectral values and a v k . However, when tp* applied to invariant elements with 
respect to the above Weyl group, the ambiguity vanishes. 

Note also, that Schubert cocycle a w £ H* (J-"a(?i)) is invariant with respect 
to W(D) if and only if w is the shortest representative in its left coset modulo 
W(D). Such cocycles form the canonical basis of cohomology H*(T a {'H))- 

3.3.2. Schubert polynomials. To calculate the coefficients c^(a) we have to return 
back to the Schubert cocycles o~ w and express them via the canonical generators 
Xi. This can be accomplished by the divided difference operators 

di l f(Xi,X 2 ,...,X n ) H-> (41) 

Xi 2^+1 

as follows. Write a permutation w S S n as a product of the minimal number of 
transpositions Sj = (i, i + 1) 

w = s n s i2 ■ ■ ■ s ie . (42) 

The number of factors £(w) = #{z < j \ w(i) > w(j)} is called the length of the 
permutation w. The product 

is independent of the reduced decomposition and in terms of these operators the 
Schubert cocycle a w is given by the equation 

o~w — d w -i Wo (x-^ X2 ■••x n —x), (43) 

where wq — (n, n — 1, . . . , 2, 1) is the unique permutation of the maximal length. 

The right hand side of equation (|43| makes sense for independent variables xi 
and in this setting it is called Schubert polynomial S w (xi, X2, . . . ,x n ), deg S w = 
£(w). They where first introduced by Lascoux and Schiitzenberger [2T1I22] who 
studied them in a long series of papers. See |24j for further references and a 
concise exposition of the theory. We borrow from [21] the following table, in 
which x, y, z stand for x\, X2, X3. 



w 


S w 


w 




w 


S w 


w 


S w 


3210 


3 'A 

x y z 


2301 


•2 '2 

x y 


2031 


'A . '2 

x y + x z 


1203 


xy 


2310 


'2 '2 

x y z 


3021 


x s y + x 1 z 


2103 


x z y 


2013 


x 1 


3120 


x A yz 


3102 


x*y 


3012 


x J 


0132 


x + y + z 


3201 


*V 


1230 


xyz 


0231 


xy -t- yz + zx 


0213 


x + y 


1320 


x yz + xy z 


0321 


•2 , 2 , 2 

x y + x z + xy 


0312 


x z + xy + y z 


1023 


X 


2130 


x-'yz 


1302 


'2 1 '2 

x y + xy 


1032 


x 2 + xy + xz 


0123 


1 



Extra variables x n +i, Xn+2t ■ ■ ■ being added to (|43[) leave Schubert polyno- 
mials unaltered. By that reason they are usually treated as polynomials in an 
infinite ordered alphabet x = (x\, X2, ■ ■ ■)■ With this understanding every homo- 
geneous polynomial can be decomposed into Schubert components as follows 



f( x ) = ^2 ®wf ■ S w (x). 

£(w)=dc e (f) 



18 



Murat Altunbulak and Alexander Klyachko 



Applying this to the polynomial 



e{v)=i(w) 



and using Proposition |4] we finally arrive at the following result. 

Theorem 5. For the v-representability problem the coefficients of the decompo- 
sition (p*(a w ) = J2 V c w( a ) a v are gi ven by the formula 



where the tableau T is derived from equation a v k —aT 7 and the operator d v acts 
on the variables Xi , replacing x v k via specialization x v k i— ► Xt = X^ieT x i • ^ 

Take notice that this equation is independent of an ordering of the unresolved 
spectral values a v k . Indeed, Schubert polynomial S w {x v ) is symmetric in the 
respective variables a?j£, provided that w is the minimal representative in its left 
coset modulo centralizer of the spectrum a" in the symmetric group. Only such 
permutations correspond to Schubert cocycles a w £ H*(J r a ^(H 1 ^)), cf. RemarkHJ 

4. Beyond the basic constraints 

Here we use the above results to derive some general inequalities for the pure 
i/-representability problem beyond the Pauli constraint A -< v. We start with a 
complete solution of the problem for two-row diagrams, and then turn to the 
initial iV-representability problem that appears to be the most difficult one. 

4-1. Two-row diagrams. For two-row diagram v = [a,/3] the majorization in- 
equality A -< v just tells that Ai < a. As we know, for j3 ^ 1 it shapes the whole 
moment polytope, see Remark[3]to Theorem^ Here we elucidate the remaining 
case v = [N — 1, 1], and thus solve the pure ^-representability problem for all 
two-row diagrams. The result can not be extended to three-row diagrams, nor 
even to three fermion systems, where the number of independent inequalities 
increases with the rank, see Corollary [3] below. For convenience and a future 
reference we collect in the next theorem all known facts. 

Theorem 6. For a system 1~L v r of rank r > 3 with two-row diagram v = [a,/3], 
a + (3 — N all constraints on the occupation numbers of a pure state are given 
by the following conditions 

1. Basic inequality \\ < a for f3 =/= 1. 

2. Inequality Ai - A 2 < N - 2 for v = [N - 1, 1], N > 3. 

3. Inequalities Ai — A2 < 1, A2 — A3 < 1 for v = [2, 1]. 
4- Even degeneracy A2i-i = A2i for v = [1, 1]. 

Proof. We have already addressed the cases 1 and 4 in Remark [3] and Introduc- 
tion respectively. 

Case 2: Necessity. To prove the inequality Ai — A2 < N — 2 we have to put it 
into the form of Theorem [2] 



c w( a ) — d v S w (x ) |x£i->'. 



(44) 




(45) 
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This suggests the test spectrum a — (1, 0, 0, . . . , 0, — 1) and the shortest per- 
mutation v that transforms it into (1, — 1, 0, 0, . . . , 0), which is the cyclic one 
v = (2, 3, 4, ... , r). Thus we get the left hand side of the inequality. To interpret 
its right hand side N — 2, notice that the spectrum a" starts with the terms 

a v = (N - 1 , N - 1, . . . , N - f , N - 2, . . .), 

" V ' 

r-2 

corresponding to semi-standard tableaux T with first row of ones and the indices 
2, 3, . . . , r filling the unique place in the second row. Since for pure state [i = 
(1, 0, 0, . . . , 0), then the shortest permutation w that produces N — 2 in the 
right hand side of (143)) is also cyclic w — (1, 2, 3, . . . , r — 1). The corresponding 
Schubert polynomial is just the monomial 

Snj{^X ) X^X^ ' ' ' X T 2' 

This is a special case of Grassmann permutations discussed in the next n° 14.21 
Specialization x v k i— > xt of Theorem [5] transforms it into the product 

r- 1 

P(x) = Y[[(N-l)xi+Xi]. 

i=2 

Taking the reduced decomposition v = S2S3 ■ ■ ■ s r _i we infer 

cl(a) = d v P{x) = d 2 d 3 ■ ■ ■ dr-iP(x). 

The right hand side is a constant, and the operators di do not touch x\. Hence 
we can put X\ — 0, that gives 

c v w {a) = d 2 d 3 ■ ■ ■ d r - 1 (x 2 x 3 ■ ■ -x r _i) = 1. 

Since c^(a) ^ 0, the inequality follows from Theorem [2] 

Case 2: Sufficiency. By the convexity it is enough to construct extremal states 
whose occupation numbers are vertices of the polytope cut out from the Weyl 
chamber by the inequality Ai — A2 < N — 2 and the normalization Tr A = N . The 
vertices are given first of all by the fundamental weights normalized to trace N 

uj k = (N/k,N/k,...,N/k,0,0,...,0) 

s y ^ 

k 

that generate the edges of the Weyl chamber, except for loi forbidden by the 
constraint Ai — A2 < N — 2. The latter is replaced by the intersections of 
segments [a>i,a;fe] with the hyperplane Ai — A2 = N — 2 

r fe = (N - 2 + 2/fc, 2/fc, . . . , 2/fc, 0, 0, ... , 0). 

k 

Here we tacitly assume that N > 3, since otherwise aj 2 would be also forbidden. 
The same condition ensures that the system TC% is non-exceptional for k > 2, 
hence Uk are occupation numbers of some pure states by Proposition [2j 
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To deal with the remaining vertices we invoke the Dadok-Kac construction 
n° 13.2.31 and observe that the state 



4>k 



_ 1 fc fe 



V 2<i<k 



has a disconnected support and the occupation numbers Tfc, fc > 2. Here for 
clarity we write tableau T instead of the weight vector er and skip an overall 
normalization factor. 



Case 3. Here we only briefly sketch the proof that follows a similar scheme. The 
second inequality in the form A2 — A3 < N — 2 holds for all N, but it becomes 
redundant for N > 3. It can be deduced from Theorem [5] by calculation of the 
coefficient c v w (a) for the same a and w as above, but with another permutation 
v — (1, 2)(3, 4, . . . , r). Then, keeping the notations of Case 2, we get 

c v w (a) = <9 3 <9 4 ■ ■ -dr-idiP(xi,x 2 , . . . ,X r -l) 

P(xi,X 2 , ■ ■ • ,2V- 1) - P(x 2 ,Xi, . . . ,X r -l) 
— C3C4 • • ■ r _i . 

x\ — x 2 

The operators c\, k > 3 do not affect variables x±,x 2 . Therefore we can pass 
in the fraction to the limit X\, x 2 — > equal to (N — 2)^3X4 • • • x r -i, that gives 
cl(a) = N-2^0. 

To prove sufficiency of the above inequalities we again have to look at the 
vertices of a polytope cut out of the Weyl chamber by the constraints \i — \ 2 < 1 , 
A2 — A3 < I, Tr A = 3. This time, along with Uk, k > 3 and Tfc, k > 2, there are 
vertices of another type 



T] k = (I + I /fc, I + I /fc, 1/fc, 1/fc, . . . , 1/fc, 0, 0, . . . , 0) 

v. ^ 

fe 

for fc > 3. They represent occupation numbers of the following states with dis- 
connected support 

yj k = VfcTT gp + V2jp+ W 1 - ° 

3<i<k 

Remark 5. Two-row diagrams naturally appear in description of bosonic systems, 
like photons where polarization plays role of spin. Representation with diagram 
[ip can be applied both for bosons and fermions. In this case we calculated 
all constraints on the spin and orbital occupation numbers for small ranks, see 
n° 16.11 It appears that the constraints are stable and independent of the rank. 



4-2. Grassmann inequalities. Let's return back to the initial pure iV-represent- 
ability problem for system A N 7i r and consider a constraint on its occupation 
numbers with 0/1 coefficients 



Ajj + Xi 2 H h X ip < b, 



(46) 
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called Grassmann inequality. For example, all constraints Q for system A 3 7^7 
are Grassmannian. We assume that the Grassmann inequality is essential, mean- 
ing that it defines a facet of the moment polytope. Then it should fit into the 
form of Theorem [2] with 

o = ( l,l , . ..,1 ,0,0,...,0) 
v 

and the Grassmann permutation or shuffle 

v = [ii,i 2 , . . ■ ) i p ,h,j2, ■ ■ -Jg] ■= [I, J], (47) 

where / and J are increasing sequences of lengths p and q, p + q — r. This is the 
shortest permutation that produces the left hand side of inequality (|46|) . Our 
terminology stems from the observation that for the test spectrum a the flag 
variety T a {TL) reduces to the Grassmannian Gv^iTL) consisting of all subspaces 
in TL of dimension p and codimension q. 

It is instructive to think about Grassmann permutation v = [I, J] geometri- 
cally as a path r connecting SW and NE corners of p x g rectangle, with fc-th 
unit step running to the North for k 6 I and to the East for k £ J. The path 
cuts out of the rectangle a Young diagram 7 at its NW corner. We'll refer to 
I and J as the vertical and horizontal sequences of the diagram 7 C p x q and 
denote the corresponding shuffle by v 1 — [I, J]. The length of the shuffle w 7 is 
equal to the size |7| of the diagram 7 and its Schubert polynomial reduces to 
the much better understood Schur function 

Observe that 7 P -fc+i = ik — k, and the size of the Young diagram 7 related to 
its vertical sequence by the equation 

M= E ^- k )- ( 48 ) 

l<k<p 

To get the strongest inequality (|46]) we chose w to be cyclic^ permutation 

w= (1,2,..., £+ 1) = [2,3,... ,£+ l,l,£ + 2,£ + 3..., r] 

of length £ = t(v) = \-f \ for which the right hand side b = (A N a)e+i of P5|) is 
minimal and equal to I + 1-th term of the non-increasing sequence 

A N a = {a K ■= a kl + a k . 2 -\ h a kN \ 1 < k± < k 2 < ■ ■ ■ < k N < r} 1 . 

The sequence consists of nonnegative numbers m each taken with multiplicity 

(')( - V 

\m ) \N — m J 

Recall that w also should be the minimal representative in its left coset modulo 
stabilizer of /\ N a. For the cyclic permutation this amounts to the inequality 



e Actually w is always cyclic for an essential pure y-reprcscntability inequality. We'll address 
this issue elsewhere. 
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(A N a)e > (A N a)e + i = b, which tells that the first I terms of A N a contain all 
the components bigger than b. The number of such terms is bounded by the 
inequality 

m>b v 7 v y 

To avoid sporadic constraints, assume that the inequality we are looking for is 
stable, i.e. remains valid for arbitrary big rank r. Then the left hand side should 
be linear in q = r — p and the sum contains at most two terms: m = N and 
m = N — 1. Thus we end up with two possibilities 

1. b = N — 2, p — N ~ 1, l — r— p, that gives the inequality 

A n + \ t2 + ■ ■ ■ + \ lN _ x <N-2, (50) 

witn Sfc(*fc -k) = r-p. 

2. b = N - 1, p > N, £ = (£), that gives the inequality 

A n + A t2 + • • • + A tp <N-1, (51) 

with £*(**-*) = (&)■ 

We will refer to them as the Grassmann inequalities of the first and second kind 
respectively. For the inequalities of the first kind the sum ^2 k (ik — k) = r — p 
increases with the rank, and therefore some of the involved occupation numbers 
should move away from the head of the spectrum. In contrast, the constraints 
of the second kind deal only with a few leading occupation numbers that are 
independent of the rank. We analyze them below for p — N + 1 and postpone 
a more peculiar first kind to the next section. The final result is that these 
inequalities actually hold true with very few exceptions. 

The cyclic permutation w is a special type of shuffle with column Young 
diagram of height I. The corresponding Schur function is just the monomial 

S w (y) = 2/12/2 ■■■yt- 

Applying to S w the specialization of Theorem [5] we arrive at the product 
P{x) = Yl {xk x +x k2 ^ hfffcjy) = y^c 7 5 , 7 (xi,x 2 , ■ ■ ■ ,x p ). (52) 

l<fei<fe 2 <---<fcjv<J9 7 

Being symmetric, it can be expressed via Schur functions and, by Theorem [2l 
each time Sy(x) enters into the decomposition with nonzero coefficient c 7 =/= 
we get inequality 

Xn+K + '-' + K <N-1, (53) 
where i\ < ii < ■ ■ ■ < i p is the vertical sequence of Young diagram 7 C p x q, 

M = (&)• 

The product P(x) represents the top Chern class of the exterior power A N £ P 
of the tautological bundle £ p on Grassmannian Gr* and the decomposition (|52| 
has been discussed in this context [20]. However, known results are very limited. 
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Example 3. For N = 2 and any p > N the product 

P{x) = J][ (xj + Sj) = S s (x 1 ,x 2 , ■ ■ ■ ,x p ) 

l<z<_7<p 

is just Schur function with triangular Young diagram 5 = \p — l,p — 2, . . . , 0], 
see |25j . This gives for two fermion system A 2 H the inequality 

Ai + A 3 + A 5 + A 7 • • • < 1, (54) 

that, due to the normalization J^. A; = 2, degenerates into equality and implies 
even degeneracy A21-1 = A2; of the occupation numbers. 

On the other hand, for arbitrary TV and minimal value p — N we get 

P{x) = x\ + x 2 + ■ ■ ■ + x N = So{x)- 

The vertical sequence of the one-box diagram □ gives a nontrivial inequality 

Ai+A 2 H h Xn-i + Ajv+i < N - 1 (55) 

that forces N-th electron into N-th orbital, when the preceding orbitals are fully 
occupied. We improve it below. 

To the rest of this section we focus upon the next case p = N+ 1 that provides 
an infinite series of inequalities. Observe that in this setting a row diagram 7 of 
length N + 1 = (j^j produces a false inequality 

Ai + A 2 + --- + Ajv + A 2 iv+2<iV-l, (?) (56) 

that fails for a coherent state given by one Slater determinant ei A e 2 A . . . A ejv. 
Similarly, the column inequality 

A 2 + A 3 + . . . + Aat+2 < N - 1 (?) (57) 

fails for even N. Indeed, in this case the system A n Hn+2 C l\ N TL r is non- 
exceptional and hence, by Proposition [21 the spectrum 

X=rf^{ N,N,...,N ,0,0---,0) 

N+2 

represents legitimate occupation numbers violating the inequality. 

Quite unexpectedly, all the other diagrams produce a valid constraint. In 
plain language the result can be stated as follows. 

Theorem 7. The occupation numbers of N -fermion system f\ N 7i in a pure state 
satisfy the following constraint 

K + X i2 + ■ ■ ■ + X iN+1 <N-1 

each time ^2 k (ik — k) — N + 1, except for inequality (|56p and, for even N, 
inequality (|57p. 
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Proof. For p — N + 1 the decomposition (|52l) takes the form 

P(x) = Yl {xi+x 2 -\ hxH \-x N+1 )= Yl {°i-Xi) 

l<i<N+l \<i<N+l 

= £ (-l) fc crf +1 ~ fc cr fc = y^c 7 5 7 (a:i,a;2, ■ ■ .,x N+1 ), 

0<fc<JV+l 7 

where (Jk{x) — S^k^(x) are elementary symmetric functions, or what is the same 
Schur functions for the column diagram [1*]. 

For Young diagrams r C 7 denote by tfy/r) the number of standard tableaux 
of skew shape 7/r. Then 

c 7 = Et- 1 )"^/^])- ( 58 ) 

fe>0 

Indeed, the coefficient at 5 7 in crf r+1_fc (Tfc = S^ +1 ~ k S^kj is equal to the number 

of ways to build 7 from the column diagram [l fc ] by adding cells one at a time. 
Numbering the cells in the order of their appearance gives a standard tableaux 
of shape 7/[l fc ] that encodes the whole building process. Thus the coefficient is 
t(7/[l fc ]) and the equation (|58|) follows. 

For a column diagram 7 we infer from the last equation 



JV+l 



= £(-D* = C 



k=0 



N = mod 2, 
N = 1 mod 2. 



Henceforth we assume that 7 is not a column. Let's combine successive even and 
odd terms of the sum (l58l) 



]>>( 7 /[l 2 *])-t(7/[l 2l+1 ])]. 



i>0 



We claim that 



i(7/[l fe ])-^(7/[l fc+1 ])=i(7/[2 ! l fc - 1 ]) 



(59) 



(60) 



where meaningless terms understood as zeros, e.g. the right hand side for k = 0. 

Indeed, the building process can be described as an extension of the partially 
filled tableau 



to a full standard tableau of shape 7. One can put the number k + 1 either just 
below k or next to 1. For the first choice the number of ways to complete the 
tableau is t(7/[l fe+1 ]), while for another one the number is £(7/[2, l fc_1 ]). Hence 
t( 7 /[l fc ])=*(7/[l fe+1 ]) + *(7/[2,l fe - 1 ]). 
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Combining the last two equations we arrive at the following representation of 
the coefficient c 7 as a sum of nonnegative terms 



^tfr/pU 24 - 1 ]). 



(61) 



i>0 



For a row diagram all terms vanish, while otherwise t(-y/[2, 1]) ^ 0. Hence c 7 > 
if the diagram is neither a row nor a column. The result now follows from 
Theorem H □ 

Example 4- For TV = 3 the theorem gives four inequalities listed below together 
with the corresponding diagrams 



A 2 + A 3 + A 4 + A 5 < 2, 



Ai + A 2 + A 5 + A 6 < 2, 



Ai + A 3 + A 4 + A 6 < 2, 



Ai + A 2 + A 4 + A 7 < 2. 



(62) 



They are valid for arbitrary rank r and give all constraints on the occupation 
numbers for r < 7. 

Observe also an improved version of the inequality (|55p 

Ai + A 2 + • • ■ + Ajv-i + X N+1 + X 2N +i <N-l, (63) 

coming from the diagram [N, 1] , and another inequality 

A 2 + A 3 + h Aw+2 < N- 1, 

originated from a column diagram and valid only for odd N. 

Remark 6. We have considered above only Grassmann inequalities of the lowest 
levels p — N, N + 1. The higher levels provide further improvements. For ex- 
ample, the inequalities (|55p and (|63[) are just the first terms of an infinite series 
corresponding to increasing values of p 



X tl + \ t2 + A l3 + • • • + X ip < N - 1, 



(64) 



where ik = k + (jy-Jjj) ■ For N — 2 this gives the inequality (|54[) and the double 
degeneracy of the occupation numbers, while for N = 3 we get the inequality 



Ai + A 2 + A 4 + A 7 + An + A 



16 



< 2, 



where the differences between the successive indices are natural numbers 1,2,3, 4, 
The details will be given elsewhere. 
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4-- 3. Grassmann inequalities of the first kind. Formally we have such an inequal- 
ity 

K +A i2 +--- + X iN _ 1 <N-2 (65) 
each time the Schur function S~, = S v enters into the decomposition 

P(x) = Yf_ (xi+x 2 -\ \-Xjf-i = ^ c v S v (x). (66) 

N<j<r t(v)=t 

Here 7 is a Young diagram of size £ = r — N + 1 with the vertical sequence formed 
by the indices in the above inequality, and t> 7 is the corresponding shuffle. In 
contrast to the previous case, the product is not a symmetric function and its 
decomposition into Schubert polynomials is a challenge. 

Let's try a simple case of a row diagram that produces the inequality 

Ai+A 2 + --- + A A r- 2 +A r .<JV-2. (67) 

A close look shows that it fails for odd I ~r — N +1 = 2m — 1 for the spectrum 

A = C 1 ' h}l m ' l / m > ■ ■ ■ l l / m ) 

N-2 2m 

obtained by merging of the occupation numbers of the systems A n ~ 2 Hn-2 and 
A 2 7i2m, see Corollary [2] of Theorem |U Neveretheless 

Proposition 5. The inequality (|67p holds for even £ = r — N + 1. In this case 
the Schur function with a row diagram enters into the decomposition (|66j) with 
unit coefficient. 

Proof. The row diagram 7 corresponds to the cyclic permutation 

v = v 1 = (r,r - 1, ...,N,N - 1) = s r _is r _ 2 ■ ■ -Sjv-i, 

where Si = (i,i + 1) are transpositions. We have to calculate the coefficient c v 
of the decomposition (|66j) given by the equation 

c v = d v P(x) = d r -id r -2 ■ ■ ■ d N -iP(x). 

The operator d v does not affect the variables Xi, i < N — 1, so we can set them 
to zero and deal with the polynomial 

Pa{x) = JJ_ (x N -i+Xi)= ^ x^^x^x^-'-x^. 

N<i<r N<ii<i 2 <---<i k <r 

We claim that 

a i-k _ J (- 1 ) fe for i s = r - fc + s, 

@vX _ 1 x^, xi 2 • ■ • X{, < v^^J 
JV 1 10 otherwise. 

Let start with the second case i\ < r — k — £ + N — k — 1. In the following 
calculation we set to zero all variables that are not affected by the subsequent 
operators dj . With this convention we get 

— 2<9ii — 3 ' ' ' &N—1X _ \Xi 1 Xi 2 • Xi k — _-y Xi t Xi 2 • • • Xi k . (^^) 
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The resulting monomial is divisible by s^-i-invariant factor x^-iX^ that com- 
mutes with operator c^—i. Hence everything vanishes in the next step as a result 
of the action di t -\ and setting x^-i = 0. 

In the case i\ — r — k + 1 = £ + N — k the right hand side of (|69[) is just 
the product of the last k variables x r —k+i%r—k+2 • ■ ■ x r and application of the 
remaining operators dj, r — k < j < r — 1 gives (— l) fe . 

Finally, from the equation 168[) we infer 

" is even, 
is odd. 



0<A:<^ 



and the result follows from Theorem [2] □ 

Remark 7. The inequality (|67j) is most appealing for AT = 3 

Ai+A r <l, (71) 

where it supersedes the Pauli principle Ai < 1 for even r. Note that for three 
electron system one- and two-point density matrices are isospectral and therefore 
the above inequality holds for both of them. We first came across this result 
reading paper [14] , where the authors observed that if the 2-point density matrix 
of a three fermion system in state -0 S A 3 Tt r has an eigenvalue equal to one, 
then the corresponding eigenform u) £ f\ 2 TL r can't have the full rank r. This is 
trivial for odd r, since rank of to is always even. For even rank this follows from 
([7Tj) . Moreover, in the latter case the state ip £ A 3 Tt r itself has rank less than 
r. M.B. Ruskai also conjectured inequality (f7Tj) in her analysis of three fermion 
and three hole systems [33] . 

Observe of the following result, anticipated by many experts. It may appear not 
so trivial if compared with Theorems |4] and [6l 

Corollary 3. No finite set of inequalities gives all constraints on occupation 
numbers of N -fermion system A N Tt, N > 1 of arbitrary big rank. 

Proof. Indeed, a finite set Q of linear inequalities L a (X) < b a includes only 
finitely many occupation numbers Xi, i < M. Every inequality that follows 
from Q is a nonnegative combination of the inequalities from Q, the ordering 
conditions A^ — A,_i < 0, and a multiple of the normalization equation y^._ 1 A; = 
N. 

Suppose now that the inequality of Proposition [5] 

Ai + A 2 + --- + AAr_2 + A r <N-2 (72) 

can be deduced from the system Q for some r ^> M and even I = r — N + 1. 
The coefficients at Ai in the left side for i > M should come from the following 
linear combination with non-negative coefficients 

ai(A 2 - Ai) + a 2 (A 3 — A 2 ) + • • ■ + a r _i(A r — A r _i) — a r \ r — 
— \\a\ + A 2 (ai — a 2 ) + • • • + A r _i(a r _ 2 — a r _i) + A r (a r _i — a r ) 

amended with a multiple of the normalization equation. The Abel transforma- 
tion shown in the second line implies that the coefficients ai should form an 
arithmetical progression = ai + b for M < i < r, while a r = ar + b — 1 > 0. 
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Suppose now that a > 0. Then the same combination of inequalities from 
Q that produces (|72|) and the same coefficients for i < r together with a r = 
ar + b > 0, a r +\ = o(r+l) + b— 1 > would give a false inequality of rank r+1 
obtained from (|T2[) by replacing r i— > r + 1 . Recall that the inequality (f72|) /ai/s 
for odd £ = r — iV + 1. For a < a similar consideration gives a false inequality 
of rank r — 1 . □ 

Proposition [5] can be extended to two-row diagrams 7 = — fc, fe]. For three 
fcrmions this leads to the constraints 

Afe+i + A r _fe < 1, for fc + l<r-fe, (73) 

that prohibit more than one electron to occupy two complementary orbitals. 
It holds both for even and odd r for k > 0. The corresponding coefficients 
c 7 = c(£, k) of the decomposition (j6"6")l satisfy the recurrence relation c(£, k) = 
c(£ — 1, k) + c(£ — 1, — 1) and form the left half of the Pascal triangle 
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with apex at £ = —1, and 0/1 boundary condition for k — set by equation 
(|70p . We return to the Pascal recurrence relation in a more general framework 
below, see equation (f79|) . 

Observe a zero in the forth line of the Pascal triangle, corresponding to dia- 
gram In general, a column diagram should have zero coefficient, because it 
produces inequality 

Ai + A 2 + • • • + X^ e + --- + X N <N-2 (?) (74) 

that fails for a coherent state given by one Slater determinant. 

It turns out that the Grassmann inequality of the first kind (|65p holds for all 
diagrams, except for a column and an odd row. To wit 

Theorem 8. The occupation numbers of N -fermion system f\ N 7i r in a pure 
state satisfy the following constraint 

K+K+--- + Aijv-! < iV - 2 (75) 

each time X)fe(*fe — ^) = r — N + 1, except for inequality (|74p and, for odd 
I = r — N + 1, inequality ([6T)l . 

Proof. We've to show that Schur function 5 7 (x) = S Vy (x) enters into the de- 
composition 



P r (x) = (xi+.t 2 H \-xn-i +Xj) = c v S v (x), (76) 

N<j<r £(v)=e 
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provided that 7 C p x q is neither a column nor an odd row. Here p — N — 1, 
q = £ = \j\ = r — p. 

Note first of all, that the coefficients of this decomposition are nonnegative 
for v € S r and can be positive only for shuffles v — Vj. The first claim holds in 
general for the coefficients c™ (a) of Theorem [2] 

v 

since the cycle ^a 1 (<r UJ ) c T a ^Hr) is effective. Here v runs over representatives 
of minimal length in left coset modulo stabilizer of a. To include all permutations 
one has to deal with a small perturbation a that resolves multiple entries 
of a. However, since l -p^ l {u w ) C Ta^Hr) is pull back of t^ 1 (cr K1 ) C T a (^H r ) via 
natural projection tt : Ta(Jir) —* FaCHr) defined in Remark 01 then decomposi- 
tion of </ 3 a 1 ( <J iu) an d 1 ( cr ™) involve the same Schubert cycles a v . This prove 
the second claim. Let's add as a warning, that the decomposition (|76|) actually 
contains Schubert polynomials S v with permutations v (jt S r . 

The rest of the proof is purely algebraic. We'll proceed by induction on r 
keeping N fixed. For the first meaningful case r = N + 1, i = 2, as we know, 
only row diagram | | | appears in the decomposition. 

Suppose now the induction hypothesis holds for P r (x), and consider the next 
polynomial 

P r +i(x) = (xi + x 2 H h x N -i + x r+ i)P r (x) 

= {xi+x 2 + h xn-i + x r +i) c v S v (x). (77) 

£(y)=t 

We can find its Schubert components using a version of Monk's formula 
(atxi + a 2 x 2 + a 3 x 3 ■ ■ ■ )S v (x) = ^ i a i ~ a i) S vt tj , 

where ty = i < j < 00 is a transposition, see [241 p. 86]. For a typical 

term of ([77]) this gives 

(x\ + x 2 + ■ ■ ■ +xn-i + x r+ i)S v 

E S ««- E sg*(r + l-j)S vtj , T+1 , (78) 

\<i<N<j^r+l N<j^r+l 

where the sums include only those transpositions t for which i(vt) — £(v) + 1. 
We are interested in the terms u 7 — vt € S r +i that are shuffles coming from a 
Young diagram 7 C p x (£ + 1) of size ^ + 1. Let's single out the row diagram 
for which Proposition [5] gives the coefficient c 7 . The remaining shuffles it 7 do 
not move the last index r + 1, and therefore permutation v — u^tij has a bigger 
length than u 7 for j > r+1. Hence a non-row Schur component S 7 in (|78[) comes 
from the sum 

E ^ vt 'i 
\<i<N<j<r 

for v = u-ytij, £{v) = £{u~ t ) — 1 = — 1. Then u G S r , and S v (x) enters into 
decomposition (|76[) only for a shuffle i; = u T . In this case the relation v T = u^Uj 
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just means that r is obtained from 7 by removing a cell. As a result, we arrive 
at the recurrence relation 

c 7 = ^2 c r , (79) 

7 / T— cell 

that holds for all non-row diagrams 7. This implies that c 7 > if one can obtain 
an even row from 7 by removing cells one at a time from a non-row diagram. 
This can be done for any diagram different from a column or an odd row. The 
inequality (175|) now follows from Theorem O □ 

Example 5. For four fermion system f\ A TL r the theorem gives inequality 

A( + X 3 + A fe < 2, 

that holds for odd rank r > 7 and pairwise distinct indices satisfying equation 
i+j + k = r+3. For even r one has to exclude the row inequality Ai +A2 + A r < 2. 

For two-row diagrams equation (|79|) amounts to the Pascal recurrence relation 
discussed in Remark [7] In general, it allows to get an explicit formula for the 
coefficient c 7 that is surprisingly similar to the one given in the proof of Theorem 
[8j where we borrow the notations. 

Corollary 4. 

c 7 = E(- 1 ) fe *ww) = E^/t 2 *' !])> ( 8 °) 

k>0 i>0 

where the second equality holds for diagrams 7 different from rows and columns. 

Proof. Applying equation the (|79p recurrently in conjunction with Proposition 
[5] we find out that c 7 is equal to the number of ways to obtain an even row from 
7 by removing cells one at a time from a non-row diagram. If 7 is not a row or 
a column, then the last step in the process will be [2i, 1] 1— * [2i]. Encoding the 
process by the standard tableaux, we arrived at the second formula. The first 
one follows from the identity t(7/[2i, 1]) = t(j/[2i}) — t(j/[2i + 1]), cf. the proof 
of Theorem [8j and holds for all diagrams. □ 

5. Connection with representation theory 

The solution of i/-representability problem suggested by Theorem [2] is not fea- 
sible, except for very small systems. For example, for four fermions A Ti-s we 
confront with an immense symmetric group of degree (^) = 70. Besides, listing 
of the extremal edges for systems of this size is all but impossible. A represen- 
tation theoretical interpretation of the ^-rcprcscntability discussed below often 
allows to mollify or circumvent these difficulties. 

Let's consider a composition of the Schur functors TL 1— > H." called a plethysm 

[H U Y= m \K x - ( 81 ) 

|a|=m-IH 

It splits into U(7i) irreducible components 7i x of multiplicity mft. It is instructive 
to treat the diagrams A and fi as spectra. We are interested in their asymptotic 
behavior for ^ and — > 00. Therefore we normalize them to a fixed size 
fi = /VH, A = X/\fi\, so that Trfi = 1 and Tr A = N = 
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Theorem 9. Every time mf^ ^ the couple (A, p,) belongs to the moment poly- 
tope of the system 7i v , i.e. there exists its mixed state p v of spectrum fx, with 
occupation numbers A. Moreover every point of the moment polytope is a convex 
combination of such spectra (A, p.) of a bounded size < M < oo. □ 

The theorem is a special case of Mumford's description of the moment polytope, 
see his appendix in [28) . It also holds in more general Berenstein-Sjamaar settings 

eh. 



5.1. Practical algorithm. For a fixed M the convex hull of the spectra (A, p) from 
Theorem [9] gives an inner approximation to the moment polytope, while any set 
of inequalities of Theorem [2] amounts to its outer approximation. This suggests 
the following approach to the mixed ^-representability problem, which combines 
both theorems. 

1. Find all irreducible components H x C [H V Y for \p\ < M. 

2. Calculate the convex hull of the corresponding spectra (A, ju) that gives an inner 
approximation V^} C V for the moment polytope V . 

3. Identify the facets of V 1 ^} that are given by the inequalities of Theorem [2] They 
cut out an outer approximation T 5 ™* D V. 

4. Increase M and continue until Vf} = ■ 

The algorithm became practical by generosity of the authors of LiE package [I], 
who made it publicly available. It allows to handle plethysms efficiently. We also 
benefit from Convex package by Franz [13j , who apply a similar approach to the 
quantum marginal problem for three qutrits [12pi7j. 

One can incorporate in the algorithm additional constraints on spectrum of 
the mixed state p v . In many problems this is just a restriction on the rank 
rk p v < p, that bounds the number of rows of p. For example, a pure state 
p v = \tp)(ip\ has rank one, the corresponding diagram p = [m] reduces to a row, 
and the plethysm amounts to the symmetric power S m (TL v ). More generally, for 
spin-orbital occupation numbers of a system of electrons of total spin J, we have 
to deal with mixed states of rank 2 J + 1, see Corollary [1] to Theorem [5J and 
respectively with the diagrams p of at most that height. 



5.2. Particle-hole duality. Here is another application of Theorem[9] Recall, that 
we arrived at the i/-representability problem from the spin-orbital decomposi- 
tions (|17[) of n° [3J In this setting the Young diagram v comes together with a 
rectangular frame r x s D i>, where r and s are dimensions of the orbital and 
spin spaces respectively. Let v* be the complementary diagram to v in the frame 
r x s, that is v* = s — v r +i-i- One can think about the representation as 
describing the holes of the system TC". These are dual systems with a natural 
pairing Ji" ® Ji v r — » Ti.^ xs = det(H r )® s , that can be extended to a pairing of 
the plethysms [H u r f <8> [UfY ~> det(H r )® sm , where m = \p\. The latter duality 
means that if Ti x is a component of then TL X is a component of \H v r ] M 

of the same multiplicity. Here A* is the complementary diagram to A C r x sm. 
In view of Theorem [9] this implies 
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Corollary 5. The moment polytope of the hole system 7i" is obtained from 
the moment polytope ofH v r by the transformation (A, /i) i— > (A*,/z), where A* = 
s — A r+ i_i. □ 

6. Analysis of some small systems 

Here we take the challenge to explore all the constraints on the occupation 
numbers. This is clearly a mission impossible. It moves us from a garden of the 
carefully selected species we dealt with in the preceding sections, into the midst 
of a wild jungle with no order or end in sight. 

To succeed in this environment we try the algorithm n° 15.11 first. However, 
due to computer limitation, it can be accomplished only for very small systems. 
For the pure 7V-representability problem these arc the systems for which Borland 
and Dennis made their prophesy 35 yeas ago [3] . To move further we use any tool 
available, from a clever guess to a numerical optimization. The final outcome of 
this endeavour are all the constraints for the systems of rank not exceeding 10. 
For r < 8 we provide a rigorous proof below. We also have a proof for system 
A 3 Hg based on other ideas, not discussed here. For the remaining cases the 
constraints are complete only beyond a reasonable doubt. To resolve the doubt 
one has to verify independently that the vertices of the constructed polytope are 
legitimate occupation numbers. We did this using a variety of methods for most 
of the vertices, but some still evaded all the efforts. For the latter we resort to 
the numerical optimization to check that they indeed can be approached very 
closely within the moment polytope. The biggest system we treated A 5 7iio is 
bounded by 161 inequalities. 

We are ready to bet a bottle of decent wine for every additional essential 
constraint found. 



6.1. Spin and orbital occupation numbers. Let's start with a simple example of 
constraints on spin fj, and orbital A occupation numbers for a system of three 
electrons of the total spin J = 1/2. By Corollary [T] to Theorem [2] the problem 
is equivalent to mixed z^-representability for v = |ip and Spec p v = (/ii,/i 2 ). A 
calculation based on the algorithm n° 15.11 shows that the constraints amounts 
to 5 inequalities 

Ai-A 2 <l + ^ 2 , A 2 -A 3 <l + /i 2 , Ai-A 3 <2-^ 2 
A,-A 2 -A 3 <1, 2Ai - A 2 + A 4 < 4- /i 2 , 

that apparently are independent of the rank. We test them for r = 4, 5. Recall 
that A and fi are arranged in the non-increasing order and are normalized to the 
traces 3 and 1 respectively. 

6.2. Pure N -representability. The known solution for two fermions, together 
with the particle-hole duality n° 15.21 bound the pure N— representability prob- 
lem to the range 3 < TV < r/2. For rank r < 8 this leaves us with systems /\ 3 H.e, 
/\ 3 H 7 , A 3 H 8 , and A 4 7i 8 . 

For three of them A 3 7Y 6 , A 3 H.7 and A 4 7Y 8 the algorithm n° EH] runs flaw- 
lessly and terminates at M = 4, 8, 10, respectively. The independent constraints 
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Inequalities 


v e 5 6 


W 6 520 


CM 


Ai + A 6 < 1 


(2 6 5 4 3) 




1 


A 2 + A 5 < 1 


(1 2 5 4 3) 


(1 2 3 4 5) 


1 


A 3 + A 4 < 1 


(1 3)(2 4) 




1 


A4 < A 5 + Ag 


(1 4 3 2) 


(1 2 3 4) 


1 



Table 1. JV-rcprescntability inequalities for system A 



Inequalities 




v 6 SV 


w e S 35 




A 2 + A 3 


+ A 4 


+ A 5 < 2 


(1 2 3 4 5) 




1 


Ai + A 3 


+ A 4 


4- A 6 < 2 


(2 346 5) 


(1 2 3 4 5) 


1 


Ai + A2 


+ A 4 


4- A 7 < 2 


(3 4 76 5) 




1 


Ai + A2 


+ A 5 


4- A 6 < 2 


(3 5)(4 6) 




1 


Table 2. 


A^-representability inequalities for system A s Ti.7- 


Inequalities 




v£Sg 


w G Sto 


cl(a) 


Ai < 1 






(1) 


(1) 


1 


A5 — A6 - 


A7 - 


A 8 < 


(1 5 4 3 2) 




1 


Ai — A2 — 


A 7 - 


Ag < 


(2 3 4 5 6) 




1 


Ai - A 3 - 


A 6 - 


A 8 < 


(3 45 76) 




1 


Ai — A4 — 


A 6 - 


A 7 < 


(4 5 8 76) 


(1 2 3 4 5) 


1 


Ai — A4 — 


As - 


Ag < 


(4 6) (5 7) 




1 


A3 — A4 — 


A7- 


Ag < 


(1 3 2)(4 5 6) 




1 


A2 — A4 — 


A 6 - 


Ag < 


(1 2)(4 5 7 6) 




1 


A 2 + A 3 + 


A 5 - 


Ag < 2 


(1 2 3 5 4) 




1 


Ai + A 3 + 


X 6 - 


Ag < 2 


(2 3 6 54) 




1 


Ai + A 2 + 


A 7 - 


Ag < 2 


(3 76 54) 




1 


Ai + A 2 + 


A3 - 


A 4 < 2 


(4 5 6 78) 


(1 2 3 4 5) 


1 


Ai + A 4 + 


As - 


Ag < 2 


(2 4) (3 5) 




1 


Ai + A 2 + 


As - 


A 6 < 2 


(3 5 4)(6 7 8) 




1 


Ai + A 3 + 


As - 


A 7 < 2 


(2 3 5 4)(7 8) 




1 



Table 3. 7V-representability inequalities for system A 4 Hs- 



grouped by the test spectra a, together with the coefficients c v w (a), and cycle 
decomposition of the permutations v, w are given in Tables UrO 
The remaining system A 3 7is is much harder to resolve. 

6.2.1. System /\ 3 Hs- We managed to decompose plethysm S m (A 3 Hs) up to de- 
gree m = 24, but still have had a discrepancy between the inner and the outer 
approximations to the moment polytope. Actually all facets of T^"? except for 
one, fit Theorem [21 For the remaining facet 

Ai + A 5 + A 6 > 1 (?) 

we use a numerical minimization of the linear form L(X) = Ai + A5 + A6 over all 
particle density matrices. It turns out that the form attains its minimum, equal 
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to |g, at the vertex 

^(15,15,15,15,6,6,6,6). (82) 

Adding this vertex gives a polytope V whose all facets are covered by Theorem 
[3J Thus V is the genuine moment polytope for A 3 7is given by 31 independent 
inequalities listed in Table [4] 



Inequalities v E S$ w 6 S56 c w( a ) 



A 2 + 


A3 + A 4 + A 5 < 2 


(1 2 3 4 5) 






1 


Ai + 


A 2 + 


A 4 + A 7 < 2 


(3 4 76 5) 


(1 2 3 4 5) 


1 


Ai + 


A 3 + 


A 4 + A 6 < 2 


(2 346 5) 






1 


Ai + 


A 2 + 


A 5 + A 6 < 2 


(3 5) (4 6) 






1 


Ai + 


A 2 - 


A3 < 1 


(3 4 5 6 78) 






1 


A 2 + 


A 5 - 


A 7 < 1 


(1 2 5 4 3)(7 8) 






1 


Ai + 


A 6 - 


A 7 < 1 


(2 6 5 4 3)(7 8) 


(1 2 3 4 5 6) 


1 


A2 + 


A 4 - 


A 6 < 1 


(1 2 4 3)(6 7 8) 






1 


Ai + 


A 4 - 


A 5 < 1 


(2 4 3) (5 6 7 8) 






1 


A 3 + 


A 4 - 


A 7 < 1 


(1 3) (2 4) (7 8) 






1 


Ai + 


A 8 < 1 


(2 8 76 5 4 3) 


(1 2 3 4 5 6 7) 


1 


A 2 - 


A3- 


A 6 - A 7 < 


(1 2) (3 4 5 8 76) 






1 


A 4 - 


A 5 - 


A 6 - A 7 < 


(1 4 3 2) (5 8 7 6) 


(1 2 3 4 5 6 7) 


1 


Ai - 


A 3 - 


A 5 - A 7 < 


(3 4 6) (5 8 7) 






1 


A 2 + 


A3 + 


2A 4 — A5 — A 7 + A 8 < 2 


(1 4 8 7 5) 






1 


Ai + 


A3 + 


2A 4 — A5 — Ag + A 8 < 2 


(1 4 8 6 7 5 2) 


(12 3. 


. 10 11) 


1 


Ai + 


2A 2 - 


- A 3 + A 4 - A 5 + A 8 < 2 


(1 2) (3 4 8 5 6 7) 






1 


Ai + 


2A 2 - 


- A 3 + A 5 - A 6 + A 8 < 2 


(1 2) (3 5 4 8 6 7) 






1 


Ai + 


A 2 — 


2A 3 - A 4 - A 5 < 


(3 64 75 8) 


(12 3. 


. 11 12) 


1 


Ai - 


A 2 - 


A 3 + A 6 - 2A 7 < 


(2 6)(3 4 5 8 7) 






1 


Ai - 


A 3 - 


A 4 - A 5 + A 8 < 


(2 8 5 746 3) 


(12 3. 


. 12 13) 


1 


Ai - 


A 2 - 


A 3 - A 7 + A 8 < 


(2 8 734 5 6) 






1 


2Ai 


- A 2 - 


f A 4 - 2A 5 - A 6 + A 8 < 1 


(2 43 8 5 76) 






1 


A3 + 


2A 4 - 


- 2A 5 - A 6 - A 7 + A 8 < 1 


(1 4) (2 3 8 5) 






1 


2Ai 


- A 2 - 


- A 4 + A 6 - 2A 7 + A 8 < 1 


(2 6) (3 8 7 4) 


(12 3. 


. 12 13) 


1 


2Ai - 


f A 2 - 


- 2A 3 - A 4 - A 6 + A 8 < 1 


(3 8)(4 5 7 6) 






1 


Ai + 


2A 2 - 


- 2A 3 - A 5 - A 6 + A 8 < 1 


(1 2)(3 8)(5 7 6) 






1 


2Ai 


- 2A 2 


- A 3 - A 4 + A 6 - 3A 7 + A 8 < 


(2 64 5 3 8 7) 






1 


-Ai + A 3 + 


2A 4 - 3A 5 - 2A 6 - A 7 + A 8 < 


(1 4 2 3 8 5) (6 7) 


(12 3. 


. 14 15) 


1 


2Ai - 


f A 2 - 


- 3A 3 - 2A 4 - A 5 - A 6 + A 8 < 


(3 8) (4 7) 






1 


Ai + 


2A 2 - 


- 3A 3 - A 4 - 2A 5 - A 6 + A 8 < 


(1 2) (3 8) (4 7 5) 






1 



Table 4. 7V-representability inequalities for system A 3 "H 8 . 



We are actually unhappy with employment of the numerical optimization, 
that can produce no rigorous result. Nevertheless, it provides a helpful hint 
about missed vertices. After some guesses and trials we found the state 

ip = 2 [123] + \/l0[145] + \/5[347] + V2[356] + V2[258] + 2[368] + [178], 

whose occupation numbers give the vertex (|82p. This provides a rigorous proof 
of the completeness the above constraints. Here [ijk] — e; A ej A is the Slater 
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Extremal states Vertices 



[1234] 


(1 


1 


1 


1 











0) 


[1234] + [1256] + [3456] 


(1 


1 


1 


1 


1 


1 





0) 


[1234] + [1256] 


(2 


2 


1 


1 


1 


1 





0) 


[1234] + [1256] + [1357] + [1467] + [2367] + [2457] + [3456] 


(1 


1 


1 


1 


1 


1 


1 


0) 


[1234] + [1256] + [1357] + [1467] 


(2 


1 


1 


1 


1 


1 


1 


o) 


%/2[1234] + [1256] + [1357] + [2367] 


(2 


2 


2 


1 


1 


1 


1 


0) 


>/2[1234] + [1256] + [1357] + [2457] + [3456] 


(2 


2 


2 


2 


2 


1 


1 


0) 


V / 3[1234] + \/2[1256] + [1357] + [2457] 


(3 


3 


2 


2 


2 


1 


1 


0) 


\/2[1234] + ^[1256] + [1357] + [1467] + [2367] + [2457] 


(3 


3 


2 


2 


2 


2 


2 


0) 


%/2[1234] + [1256] + [1357] 


(4 


3 


3 


2 


2 


1 


1 


o) 


[1234] + [5678] 


(1 


1 


1 


1 


1 


1 


1 


1) 


%/2[1234] + [1256] + [1278] + [1357] + [1368] 


(3 


2 


2 


1 


1 


1 


1 


1) 


[1234] + [1256] + [1278] 


(3 


3 


1 


1 


1 


1 


1 


1) 


\/3[1234] + [1256] + [1357] + [1458] + [2358] + [2457] + [3456] 


(3 


3 


3 


3 


3 


1 


1 


1) 


\/2[1234] + \/2[1256] + [1357] + [1368] + [1458] + [1467] 


(4 


2 


2 


2 


2 


2 


1 


1) 


2[1234] + \/2[1256] + [1357] + [1458] + [2358] + [2457] 


(4 


4 


3 


3 


3 


1 


1 


1) 


2 [1234] + v^2[1256] + [1357] + [1368] + [2358] + [2367] 


(4 


4 


4 


2 


2 


2 


1 


1) 


V2[1234] + [1256] + [1357] + [1458] 


(5 


3 


3 


3 


3 


1 


1 


1) 


%/3[1234] + [1256] + [1357] + [2358] 


(5 


5 


5 


3 


3 


1 


1 


1) 


[1234] + [1256] + [1278] + [1357] + [1368] + [1458] + [1467] 


(7 


3 


3 


3 


3 


3 


3 


3) 


%/3[1234] + V2[1256] + [1357] + [1368] 


(7 


5 


5 


3 


3 


3 


1 


1) 


%/3[1234] + [1256] + [1278] + [1357] + [1368] + [2358] + [2367] 


(7 


7 


7 


3 


3 


3 


3 


3) 



Table 5. Vertices of the moment polytope of A 4 Hs and the corresponding extremal states. 



determinant or, in our general notations, weight vector corresponding to the 
semi-standard tableau T transpose to [ijk]. Six triplets [ijk] in the support of tp, 
excluding one [356] typesetted in bold face, form a disconnected set. They are 
remnants of our failed attempt to produce the missed vertex by the Dadok-Kac 
construction n° 13.2.31 Extra tableau [356] in the support increases the number 
of adjustable parameters, but makes the problem nonlinear. Don't ask how the 
coefficients were found. 

For those people who don't trust a computer assisted proof we give an ex- 
tremal state for every vertex of the moment polytope for the systems /\ 3 1i.r, 
A 3 Tis, and A 4 TCs listed in Tables [5][6] They are sufficient for a computer inde- 
pendent proof, provided that one takes for granted the values of the coefficients 
cl(a) in Tables HS1 

6.2.2. Systems of rank 9 and 10. The results here are less definite. Only for 
smallest system A 3 7ig we have a rigorous justification of completeness for the 
system of 52 independent inequalities. For the next one A 4 7ig we found 60 con- 
straints, that give a polytope with 103 vertices. For all of them, except for two 

[16, 16, 16, 6, 6, 6, 6, 6, 6]/21, [20, 14, 14, 14, 14, 4, 4, 4, 4]/23, 

we have proved rigorously that they belong to the moment polytope. The re- 
maining two vertices were checked only numerically. It turns out that the same 
two vertices would provide the completeness of 125 constraints for A 4 Hio- The 
occupation numbers of the remaining systems A 3 7iio an d A 5 7iio are bounded 
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Extremal states Vertices 



[123] 


(1 


1 


1 











0:0) 


M OQl 1 ri A £l 


(2 


1 


1 


1 


1 





: 0) 


fl 93 4- fl 4^1 -L T9ilfil 4- T^^fil 




i 






j 


1 


n - f\\ 

u . V) 


^/9^1 9^1 4- M d%] 4- T94fil 


f3 




2 


o 




I 


: 0) 


[123] + [145] + [167] + [246] + [257] + [347] + [356] 


(1 


1 


1 


1 


1 


1 


1:0) 


^[123] + [167] + [246] + [257] + [145] 


(2 


2 


1 


1 


1 


1 


1:0) 


^[123] + ^[145] + [246] + [257] + [347] + [356] 


(2 


2 


2 


2 


2 


1 


1:0) 


[123] + [145] + [167] 


(3 


1 


1 


1 


1 


1 


1:0) 


\/2[123] + [145] + [246] + [347] 

V L J'L J'L J'L J 


(3 


3 


3 


3 


1 


1 


1:0) 


\/3[123] + v / 2[145l + [246] + [257] 

* L J 1 v L J'L J 1 L J 


(5 


5 


3 


3 


3 


1 


1:0) 


[178] + [368] + [258] + [567] + [347] + [246] + [145] + [123] 


(1 


1 


1 


1 


1 


1 


1:1) 


^[178] + [368] + [567] + [246] + ^[145] + V2[123] 


(2 


1 


1 


1 


1 


1 


1:1) 


^[178] + [258] + [567] + ^[246] + [145] + ^[123] 


(2 


2 


1 


1 


1 


1 


1 : 1) 


^[123] + ^[145] + [246] + ^[347] + [356] + \/2[258] 


(3 


3 


3 


3 


3 


1 


1:1) 


^[178] + v/2[567] + [347] + [246] + 2[145] + \/5[123] 


(4 


2 


2 


2 


2 


2 


1:1) 


[178] + [246] + [145] + y/2 [123] 


(4 


3 


2 


2 


1 


1 


1:1) 


[178] + [258] + [246] + [145] + \/2[123] 


(4 


4 


2 


2 


2 


2 


1:1) 


[258] + [567] + [145] + \/3[123] 


(4 


4 


3 


3 


1 


1 


1:1) 


\/2[145] + [246] + [347] + [356] + \/2 [368] 


(4 


4 


4 


4 


2 


1 


1:1] 


%/2[178] + [246] + [145] + \/2[123] 


(■ r) 


3 


2 


2 


2 


2 


1:1) 


[368] + [347] + \/2[145] + ^[123] 


(■ r) 


5 


3 


3 


2 


1 


1:1) 


2 [123] + v/IH[145] + v^[347] + v^[356] + ^[258] + 2 [368] + [178] 


(5 


5 


5 


5 


2 


2 


2:2) 


[178] + [567] + \/2[145] + %/3[123] 


(6 


3 


3 


3 


2 


2 


1:1) 


2 [123] + %/2 [246] + %/3[356] + ^[567] + 2[258] 


(6 


5 


5 


5 


2 


2 


1:1) 


\/2[178] + [258] + %/2[246] + [145] + \/§[123] 


(6 


6 


3 


3 


3 


2 


2:2) 


2-\/2[145] + ^[246] + ^[347] + ^[356] + ^[368] 


(6 


6 


4 


4 


4 


1 


1:1) 


2 v / 3[123] + v / 6[145] + v / 2[356]+2[567] + v / 3[258] + \/3[178] 


(7 


5 


5 


5 


2 


2 


2:2) 


V 2[145J + 2 [246] + [347] + [356J + V 2 [368] 


(V 


7 


4 


4 


4 


2 


1:1) 


^[246] + v^2[347] + ^[258] +2[368] +2^(178] + [124] 


(9 


5 


5 


5 


3 


3 


3:3) 


^[258] + [567] + ^[347] + ^[246] + 2 [123] 


(9 


6 


4 


4 


4 


3 


3:3) 


3 [145] + ^[246] + 3[347]+2[356] + v / 3[258] + ^[368] 


(9 


8 


8 


8 


3 


3 


3:3) 


\/2[178] + [258] + ^[246] + %/5[145] + \/5[123] 


(9 


9 


5 


5 


3 


3 


3:2) 


2 [123] + \/2~[246] + %/2~[356] + v%[567] + v%[258] + [368] 


(9 


9 


9 


9 


4 


4 


2:2) 


2-v/2[145] + ^[246] + ^[347] + ^[356] + ^[258] +3[368] 


(10:10:10:10:4:4:3:3) 


V%[178] + [347] + V2 [246] + V2 [145] + 2 [123] 


(ll:f 


j:6:5:5:5:2:2) 


\/3[178] + [258] + 2 [246] + %/2 [145] + V^[123] 


(11: 


11 


6 


6:4:4:3:3) 


\/3[178] + \/2[567] + [246] + 2 [145] +^[123] 


(12 :< 


5:6:5:5:5:3:3) 


[123] + >/3[145] +2 [347] + 2[356] + ^[258] + ^[368] 


(12: 


L2:7:7:4:4:4:4) 



Table 6. Vertices of the moment polytope of A 3 Hs and the corresponding extremal states. 
The first ten lines give the same data for A 3 H7- 



by 93 and 161 inequalities, but many vertices are still waiting a confirmation by 
non- numerical methods. 

The facets and vertices of the moment polytopes for all systems of rank < 10 
are available in a computer friendly format at |http://www.fen.bilk ent.edu.tr/| ~ 
murata/IM-Representability.zip. 
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